Astron. Nachr. / AN 328, No. 3/4, 292 — [2971 (2007) / DOI 10.1002/asna.200610732 



Surface magnetic field effects in local helioseismology 

H. Schunker 1 2 >*, D.C. Braun 3 , and P.S. Cally 1 

1 Centre for Stellar and Planetary Astrophysics, Monash University, Melbourne, Australia 3800 

2 Max-Planck-Institut fur Sonnensystemforschung, D-31791 Katlenburg-Lindau, Germany 

3 NorthWest Research Associates, CoRA Division, 3380 Mitchell Lane, Boulder, CO 80301, U.S.A. 



Received 2006 Oct 20, accepted 2007 Jan 2 
Published online 2007 Feb 28 



O 

(N 



Key words Sun: helioseismology - Sun: magnetic fields - sunspots 

Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic 
signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic field 
and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phase- 
sensitive holography demonstrates that they are significantly affected in a strong, moderately inclined magnetic field. 
This research indicates that the effects of the surface magnetic field are potentially very significant for local helioseismic 
analysis of active regions. 
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The potential influence of the magnetic field at the surface 
on acoustic waves is somewhat controversial, but there is 
a growing consensus that surface magnetic effects should 
be considered and included in local helioseismic analysis of 
active regions (Braun 1997; Woodard 1997; Zhao & Koso- 
vichev 2006; Braun & Birch 2006). The showerglass effect 
(Lindsey & Braun 2005a,b) is a recently suggested phe- 
nomenon, consisting of large amplitude and phase distor- 
tions of the surface wavefield due to magnetic fields in the 
photosphere. Lindsey & Braun (2005b) apply a correction 
for these effects based on the strong correlation between 
magnetic field strength and the phase perturbations. A pe- 
culiar enhancement of the phase perturbation is noted in 
the penumbra of sunspots, called the penumbral acoustic 
anomaly (Lindsey & Braun 2005a). The magnetic field in 
the penumbra is not as strong as in the umbra and this 
research attempts to quantify the phase perturbations pro- 
duced specifically by inclined magnetic fields which char- 
acterize sunspot penumbrae. 

This paper expands upon previous research by Schunker 
et al. (2005) and explores the effects that the inclined field 
within the penumbrae of two sunspots may have on acoustic 
waves originating below the solar surface. The analysis uses 
helioseismic holography (Sect. [2]). Generally, in the quiet 
Sun, the "ingression" (the deduced amplitude of incoming 
acoustic waves from a surrounding pupil) and the observed 
surface signal agree well (Lindsey & Braun 2005a,b). In 
magnetic regions, a deviation of the amplitude and phase 
of the incoming acoustic waves is indicated. Schunker et al. 
(2005) demonstrated that there is a clear cyclic variation of 
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the ingression phase with azimuthal angle within a sunspot 
penumbra, and the line-of-sight direction. It is suggested 
that the effect may be due to the process of mode conversion 
as discussed by Schunker & Cally (2006). It has also previ- 
ously been established that mode conversion is able to sim- 
ulate the observed acoustic absorption in inclined magnetic 
fields (CaUy & Crouch 2003; Cally, Crouch & Braun 2005). 
This has encouraging prospects for explaining the observed 
penumbral and inclined field dependencies. Initially, a fast 
acoustic wave propagates up to the surface where it encoun- 
ters the a = c layer (where a is the Alfven speed and c is 
the sound speed) and undergoes conversion to a fast (mag- 
netic) wave and transmission to a slow (acoustic) wave. The 
amount of energy devoted to each depends on the attack 
angle between the path of the ray and the magnetic field 
where the conversion/transmission occurs (where a = c). It 
is thought that the slow acoustic wave is what may be con- 
tributing to the observed effect of Schunker et al. (2005), 
and that presented here (see also Schunker & Cally 2006; 
Cally 2007). 

The sunspot analysis is done using full disk data from 
the Michelson Doppler Imager (MDI) (Scherrer et al. 1995) 
over a period of days as the sunspot crosses the solar disk. 
The data is tracked and Postel projected into data cubes for 
each day. This enables the penumbral magnetic field to be 
viewed from a varying line-of-sight at different heliocen- 
tric angles (one for each day of observation). The phase 
shift of the incoming acoustic waves is determined from 
the correlation of the incident acoustic wave (estimated us- 
ing holography) and the surface velocity. Vector magne- 
tograms obtained from the Imaging Vector Magnetograph 
(IVM) (Mickey et al. 1996), for both sunspots are used to 
determine the orientation and strength of the field in rela- 
tion to the phase perturbation. 
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2 Holography 

Based on observed surface signals helioseismic holography 
is able to detect subsurface wave speed variations (Lindsey 
& Braun 2000). The amplitudes of acoustic waves propa- 
gating through the interior are inferred at a focal point of 
particular depth, z, and horizontal position, r. 

For this purpose the focal plane is located at the surface, 
z = 0. The inferred contribution of the observed waves is 
calculated through the interior and back up to the surface. 
The ingression is the superposition of the amplitudes of the 
incoming waves at (r, z) and time f, 

H_(r,z,t) = I dt' I dV 

J J a< \r— r' | <b 

G-(\r-r'\,z,t-t')iP(r',t), (1) 

which is calculated here at z = 0. The Green's function 
G_ represents the sub-surface disturbance at (r, z, t) result- 
ing from a unit acoustic impulse originating at surface co- 
ordinates (r',0,t'). A Green's function derived in a wave- 
mechanical formalism that includes effects for dispersion 
(Lindsey & Braun 2000), is used here. The computation is 
confined to an annulus or 'pupil' surrounding the focal point 
r with inner and outer radius a and b respectively. Essen- 
tially, the ingression is what should result from the incom- 
ing acoustic waves propagating from the pupil to the focal 
point in the absence of perturbations. 

The ingression at the focus, r, is correlated with the sur- 
face velocity signal at r to gauge the effect of any surface 
anomaly. This "local ingression control correlation" (Lind- 
sey & Braun 2005a) is given as 

C(v) = {H^{r,v)i>*{r,v)) Av , (2) 

where ijj represents the temporal Fourier transform of the 
surface disturbance ip, v is the temporal frequency, and H- 
is the temporal Fourier transform of the ingression. The as- 
terisk denotes complex conjugation and the brackets indi- 
cate an average over a positive frequency range Av of 1 
mHz centered at a certain frequency (here at either 3, 4 or 5 
mHz). The effect of the surface perturbations is quantified 
by the phase of the local ingression correlation, 

5cp = arg[C(i/)]. (3) 

We employ observations of two mature, nearly axisym- 
metric sunspots. The observations of the sunspot in AR9026 
span a time from 2000 June 3 to June 12, and the sunspot 
in AR9057 from 2000 June 24 to 2000 July 2. This allows 
different line-of-sight observations of the sunspots as they 
rotate from approximately 60° East of the central meridian 
to about 60° West of the central meridian. These sunspots 
were selected as they do not evolve significantly over the 
duration of the observations. 

As Schunker et al. (2005) described, the full-disk MDI 
data were analyzed in 24 hour sets. For each day Postel pro- 
jections were made, centered near the sunspots. From the 
temporal Fourier transform, the ingression was computed 
and correlated to the surface Doppler signal (Eqs.[T]and[3]l. 
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Fig. 1 A representation of the geometry of 6 p . The 
'sunspot' has a radial vector k and a magnetic field vector, 
B, creating a plane. The line-of-sight vector, t, is projected 
onto this plane creating, t p , which makes an angle 9 p to k. 
This is essentially an independent parameter representing 
the line-of-sight viewing angle of the magnetic field. 



The pupil size is a = 20.7 Mm and b = 43.5 Mm for 
the inner and outer radii respectively. At a frequency of 5 
mHz, this pupil selects p-modes with spherical harmonic 
degrees {£) and radial order (n) between £ s» 450 (n = 5) 
and £ w 700 (n = 4). This pupil is also of a size to elimi- 
nate the majority of the active region signal when comput- 
ing the correlation inside the penumbra. Simple analysis of 
the acoustic dispersion relation near the surface shows that 
cot r\ w 2-kvRq /c£, where Rq is the solar radius, r\ is the 
propagation angle from vertical and c is the local sound 
speed. At 3, 4 and 5 mHz, and in the given range of £, 
cot r\ w 10, indicating a primarily vertical propagation. 

3 Ingression correlation phase 

The penumbral regions for each sunspot, inner and outer 
radii, are determined from MDI continuum images as 7.3 — 
16 Mm from the centre of the sunspot in AR9026 and 6 — 
13 Mm for the smaller sunspot in AR9057. To determine 
the dependence of the phase shift (the ingression correlation 
phase) 8<p, on the line-of-sight direction projected onto the 
plane of the local magnetic field, a parameter 6 P is defined 
(see Fig.|T}. The line of sight vector, t, is projected onto the 
plane containing the radial vector, k, and the magnetic field 
vector B. The projected line-of-sight vector, t p , thus makes 
an angle 9 p with k. 

The magnetic field vector is determined from IVM mag- 
netogram data. This provides the orientation and strength 
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Fig. 2 3 mHz ingression correlation phase (S<j>) in the 
penumbra of the sunspot in AR9057 plotted against pro- 
jected angle 9 p at different values of magnetic field incli- 
nation. The top panel (a) shows 7 < 42°, where the mean 
field strength is (B) = 1700 G, the middle panel (b) shows 
42° < 7 < 66°, where (B) = 1000 G, and the bottom 
panel (c) shows 7 > 66°, where (B) = 600 G. The hor- 
izontal dashed lines indicate the mean value of S(f> for each 
panel. The error bars indicate the standard deviation of the 
mean over bins of 20 measurements in 9 p . The solid line is 
a fit for all of the displayed data; the dotted line is a fit for 
the data from 2000 June 24-28; the dashed line is a fit for 
data from 2000 June 29-2. 



Fig. 3 3 mHz ingression correlation phase (Scf>) in the 
penumbra of the sunspot in AR9026 plotted against pro- 
jected angle 9 p in three bins of magnetic field inclinations. 
The top panel (a) shows 7 < 42°, where the mean field 
strength is (B) = 1900 G, the middle panel (b) shows 
42° < 7 < 66°, where (B) = 1400 G, and the bottom 
panel (c) shows 7 > 66°, where (B) = 600 G. The hori- 
zontal dashed lines indicate the mean value of 5(f) for each 
panel. The error bars indicate the standard deviation of the 
mean over bins of 20 measurements in 9 p . The solid line is 
a fit for all of the displayed data; the dotted line is a fit for 
the data from 2000 June 3-7; the dashed line is a fit for data 
from 2000 June 8-12. 



of the surface magnetic field in the sunspot. The IVM in- 
strument observes in the same line as the MDI (Ni 676.8 
nm), and hence at the same geometrical height in the pho- 
tosphere. Rotation and scaling were applied to the data to 
align with the MDI magnetograms. The IVM observations 
were made during a 28 minute interval on 2000 June 5, for 
the sunspot in AR9026, and 2000 June 28, for the sunspot 
in AR9057. 

It is desirable to investigate the acoustic properties of the 
magnetic field as functions of field strength and inclination. 
However, in sunspot penumbrae these quantities are not eas- 
ily separable since the magnetic field strength exhibits a 
close nearly linear relationship with the inclination from 
vertical (see Fig. 4 of Schunker et al. 2005). As in Schunker 
et al. (2005), the penumbrae are divided into three roughly 
even bins, of 7 < 42°, 42° < 7 < 66° and 7 > 66°. This 
corresponds to a mean magnetic field strength of 1900, 1400 
and 600 Gauss for the sunspot in AR9026, and 1700, 1000 
and 600 Gauss for the sunspot in AR9057, respectively. 

The behaviour of the ingression correlation phase with 
respect to the projected line-of-sight for both sunspots, and 



several frequency bandpasses, are presented in Figs. [2] to 
[7] The last figure corresponds to Fig. 5 of Schunker et 
al. (2005), which showed only the results for one sunspot 
with waves averaged over a 1 mHz frequency bandpass 
centered at 5 mHz. The results for two sunspots averaged 
over a 1 mHz frequency bandpass centered at 3, 4 and 5 
mHz demonstrates the consistency of the results for differ- 
ent spots and frequencies. 

The horizontal dashed line is the average of 5<fi over the 
entire range of 9 p . The error bars indicate the standard devi- 
ation of the mean over bins of 20 measurements in 9 p . The 
observed variation of 5(f> with the projected line-of-sight an- 
gle is consistent with the conversion of vertically propagat- 
ing acoustic waves to elliptical motion in the inclined mag- 
netic field (Schunker et al. 2005). To quantify this, the vari- 
ation of phase shift and amplitude (not shown) of the local 
ingression control correlations assuming elliptical motion, 
with the eccentricity, inclination and amplitude of the el- 
lipses as free parameters in the fits is modelled. The mod- 
els will be discussed in further detail elsewhere (see also 
Schunker 2006), but the results of the fits to 6<fi are shown 
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here to guide the eye. The solid line is a fit for all the dis- 
played data, the dotted line is a fit for data straddling the 
magnetogram (2000 June 3-7 for the sunspot in AR9026, 
and 2000 June 24-28 for the sunspot in AR9057) and the 
dashed line is the fit for the remaining data. These alternate 
fits are shown to assess the consistency of the results over 
different time periods, since only a single set of IVM obser- 
vations are used and it is assumed that the magnetic field is 
virtually static in time. The trends show similar properties 
for all time periods of the data. 

The phase of the ingression correlation is seen to vary 
across the line-of-sight angle, from P = —60° to 9 p = 
+60°, especially at stronger fields (or smaller inclinations). 
This variation progressively decreases at weaker fields. At 
stronger fields where the effect is most prominent, a slight 
dependency on frequency is observed in the two spots. At 3 
mHz, 8(f> increases by about 70° as 8 P increases from —60° 
and +60°. At 5mHz, the increase in 5(f) is about 110° as 
determined by the fits. 

4 Quiet- Sun control experiment 

The Evershed effect is a steady outflow of ~ 2 km s _1 seen 
to occur along the penumbral fibrils of sunspots (Evershed 
1909). Since this penumbral flow is roughly axisymmet- 
ric, its line-of-sight component may have a similar spatial 
dependence as the ingression correlation phase. For exam- 
ple, both quantities switch sign between the side of the 
penumbrae extending toward disk center and the side ex- 
tending away from disk center. The question naturally arises 
as to whether the variation of phase may have non mag- 
netic cause, and is perhaps related to the observed line-of- 
sight component of the Evershed flow. In fact, the Evershed 
flows increase with radial distance from the sunspot cen- 
tre. The effect presented here, actually decreases with radial 
distance, and is most predominant close to the umbra. In 
addition, we do not expect a correlation between 8(f) and the 
line-of-sight velocity, since the ingression is computed over 
a full annular pupil and should not be sensitive (to first or- 
der) to solar flows. Nevertheless, the possibility of some un- 
known physical relationship or instrumental, measurement, 
or analysis artifact should be tested. 

A control experiment was performed to assess any de- 
pendence of the ingression control correlation phases in the 
quiet Sun with the line-of-sight component of supergranu- 
lation, as determined from averages of full disk MDI Dopp- 
lergrams. This analysis must be done with care to remove 
all magnetic, including network, regions, which are charac- 
terized by downflows and may bias the result. In addition, 
large-scale spatial trends must be removed from the con- 
trol correlation phases and the Dopplergrams, since only 
possible correlations over the size of a typical penumbrae 
are relevant to interpretations of the sunspot results. Figure 
[8] shows the resulting averaged ingression control correla- 
tion phase with line-of-sight velocity over a wide area of 
the solar disk. The solid diagonal line represents the ex- 
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Fig. 4 The same as for Fig. |2l except at 4 mHz. 
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Fig. 5 The same as for Fig. [3j except at 4 mHz. 



pected phase shifts of supergranular flow if the hypothe- 
sis has merit. This assumes that the relationship between 
the phase shifts and flows is linear and must extrapolate 
to values consistent with what is observed within the two 
sunspots. The quiet-Sun phase shifts are clearly not consis- 
tent with this expectation. We believe that the Evershed ef- 
fect does not significantly contribute to the effect described 
here and in Schunker et al. (2005). 
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Fig. 6 The same as for Fig. [2] except at 5 mHz. 

5 Discussion and conclusions 

The preliminary results first presented by Schunker et al. 
(2005) deserve confirmation and elaboration due to the im- 
plications of the results. Evidence for mode conversion in 
sunspots is critical for our understanding of the MHD of 
waves and sunspots (magnetohelioseismology). In addition, 
the surface effects of the magnetic field may alter the re- 
sults of sub-surface imaging using helioseismology within 
active regions. Here, the results of Schunker et al. (2005) 
have been confirmed for a second sunspot, and it is also 
shown that the results are similar for all observed frequen- 
cies. 

A frequency dependency is supported by the theory of 
Schunker & Cally (2006). A wave of lower frequency will 
experience the upper turning point at a lower depth than a 
higher frequency wave. Therefore it will not be strongly 
affected by the magnetic field. In regions of stronger field 
strengths, corresponding to the inner penumbrae, acoustic 
waves at 5 mHz are affected more by the magnetic field 
as seen in the observations presented here. However, it is 
not entirely clear why similar trends are not observed in the 
other regions of the penumbrae. 

The effect is most prevalent in strong magnetic fields. 
Here, the observed variations correspond to travel time per- 
turbations of approximately 1 minute at all observed fre- 
quencies between 3 mHz and 5 mHz. In comparison with 
travel times used to deduce sound speeds below sunspots 
this is considerable. Zhao & Kosovichev (2006) show ev- 
idence for a similar (but smaller, 0.4 minute) variation of 
travel times with azimuthal angle around a sunspot penum- 
bra. However, their measurements were averaged over the 
entire penumbrae. 

© 2007 WILEY- VCH Verlag GmbH & Co. KGaA, Wdnhcim 





1 


I (a) 7 < 42° 























-1 
-2 


1 










1 


I (b) 42°< 7 < 


66° 






(rod) 











■ . y/ 


*o 


-1 
-2 


;• 










1 


(c) y ■ (>()' 















i _r _ 








-1 


'•■/ ' 










-2 











-100 -50 50 100 

9„ (deg) 



Fig. 7 The same as for Figj3]except at 5mHz. 

As it is shown here, the effect is highly dependent on 
the field strength, and is diminished in weaker and more 
inclined magnetic fields. However, whether this is due to 
the magnetic field strength or inclination is unclear, since in 
a sunspot the two properties are inseparable. Significantly, 
Zhao & Kosovichev (2006) also find a lack of variations in 
travel-times around the penumbra when they use MDI con- 
tinuum intensity, rather than Dopplergrams. We note that the 
interpretation of Schunker et al. (2005), namely that these 
variations are caused by elliptical motion, would predict no 
variation with line-of-sight angle of scalar quantities such 
as the continuum intensity. 

The Evershed effect may be eliminated as a major cause 
for the effect seen in the ingression control correlation. The 
ingression phase shift is larger closer to the umbra, whereas 
the Evershed effect is stronger close to the outer boundary 
of the penumbra. Line-of-sight supergranulation velocities 
do not show significant correlation with ingression phase 
shifts, which leads to the belief that what is being seen in 
penumbrae is likely a superficial variation in the ingression 
correlation phase with line-of-sight angle in magnetic fields. 

These results have important implications for helioseis- 
mic calculations within active regions. The fact that there 
is such a dependence of the phase-shift with the line-of- 
sight suggests that this is predominantly a surface effect. 
Zhao & Kosovichev (2006) have argued that time-distance 
inversions do not significantly change with the variation of 
mean line-of-sight angle of sunspots from day to day. How- 
ever, their inversions do not explicitly include or test effects 
of unresolved surface terms and so cannot directly answer 
the question of how the possible inclusion of those terms 
might change existing inferences about subsurface condi- 
tions. Theory and observations of waves in active regions 
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Fig. 8 The 5 mHz ingression phase as a function of line-of-sight velocities in the quiet sun. The diagonal solid line 
represents the expected supergranular flow. 



will certainly aid in understanding and ameliorating the ef- 
fects of surface magnetic fields with the goal to improve 
helioseismic interpretations of sunspot structure. 
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